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Abstract: A method yielding simple re- 
lationships among bilateral birth-and-death 
processes is outlined. This allows one to re- 
late birth and death rates of two processes 
in such a way that their transition probabili- 
ties, first -passage-time densities and ultimate 
crossing probabilities are mutually related by 
some product-form expressions. 



1. Introduction 

Birth-and-death processes are a powerful tool 
to describe stochastic models arising in popu- 
lation dynamics (see for instance Ricciardi [7]) 
and also in a large variety of applied fields, 
such as adaptive queueing systems and neu- 
rophysiology. In particular, their transition 
probabilities and first-passage-time densities 
play a relevant role in many applied contexts. 
Hence, obtaining closed form expressions for 
these functions is an important task. Un- 
fortunately, apart from a few cases, such ex- 
pressions are rarely encountered in the liter- 
ature. Furthermore, procedures which have 
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been successfully exploited for continuous- 
state processes are not suitable for point pro- 
cesses due to the discreteness of their state- 
space. 

The aim of this note is to give an exten- 
sion of a method leading to direct transfor- 
mations between birth-and-death processes 
(cf . Di Crescenzo [3] or [1] ) to the doubly infi- 
nite state-space case. Here we consider a pre- 
assigned birth-and-death process X t whose 
state-space is the set of integers. In Section 2 
we show how to construct a new birth-and- 
death process X t in such a way that the rates 
of the two processes are mutually related. 
As a consequence, their transition probabil- 
ities are related by a simple product-from re- 
lation. In Section 3 it is shown that also 
the first-passage-time densities of the two 
processes are similarly connected. We point 
out that the transformation from X t to X t 
can be viewed clS cL method to construct new 
stochastic models. Indeed, making use of such 
method, starting from a birth-and-death pro- 
cess with constant rates, in Section 4 a new 
process is obtained. 



2. Transition probabilities 

Let {X t ; t > 0} be a birth-and-death process 
whose state-space is the set of integers Z = 
{. . . , —1, 0, 1, . . .}. As usual, for all n G Z we 



denote by A n and fi n the birth and the death 
rates of X t , i.e. 

A n = lim - Pr{X i+5 = n + 1 1 X t = n}, 
fi n = lim ^ Pr{X t+s = n - 1 1 X t = n} 

for all n G Z. We also assume that rates 
A n and /i n are positive, so that the birth- 
and-death process has no absorbing or reflect- 
ing states. In this case X t is said to be a 
bilateral birth-and-death process (cf. Ismail 
et al. [5]). We also assume that X t is simple; 
hence, the set of rates {A„, fi n ] uniquely deter- 
mines the birth-and-death process. Accord- 
ing to Callaert and Keilson (cf. [2]), a bilateral 
birth-and-death process is simple if and only 
if the two component birth-and-death pro- 
cesses, obtained by locating at n = a bound- 
ary reflecting in both directions, are simple 
(necessary and sufficient conditions are also 
given in Pruitt [6]). 

As X t is simple, the transition probabilities 
Pk,n(t) = Pr {X t = n | X = k} are the unique 
solution of the following system: 



TrPM^) = K-lPk,n-l{t) 

at 



(1) 



- (A n + fi n ) p k ,n(t) 

+ fJ, n +iPk,n+i{t) (neZ). 
This has to be solved with initial conditions 



Pfc,n(°) 



1 if n = k 
if n ^ k, 



(2) 



where k G Z is the initial state. 

We shall now see that, under suitable as- 
sumptions, there exists a bilateral birth-and- 
death process X t having state-space Z whose 
rates are obtained from those of X t and such 
that the transition probabilities of the two 
processes are mutually related by a product- 
form relation. 



Theorem 1 Let {. . . , U-i, Uq, . . .} be a 
strictly monotonic sequence of positive num- 
bers satisfying 

V n +lK - V n (\ n + fi n ) + Vn-l^n = (3) 

for all n G Z. If a bilateral birth-and-death 
process X t has rates 



A„ = A n fi n = fi n [n G ZJ 



(4) 



then its transition probabilities are given by: 
PkA t ) = —PkA t ) (k,neZ). (5) 
Proof. Let 

Pin(t) = -PkAt) (k,neZ). 
Making use of ([1]), for all n G Z we have 

= X n -l-^P*k ,n-l(*) 
v n -\ 



- (A n + fi n ) P* k>n (t) 

From relations ([3]) and (J4]) it follows 

An ~r fl"n — r fi n 

A n -|- fi n , 

An-l — An-1 , 

Vn-1 

A*n+1 — f'n+l ■ 



V, 



n+1 



Hence, the previous system becomes: 
d 



dt 



Pk,n(t) = K-lPt,n-l(t) 

~ &n + Jln)pl,n(t) 
+ fin+lP*k, n+1 (t) (nGZ). 



These are the forward equations of a bilateral 
birth-and-death process X t , characterized by 
rates \ n ,p, n . The corresponding initial condi- 
tions are: 



lim — Pk, n (t) 

"kn "kn ■ 



Due to the uniqueness of the solution of the 
forward equations, probabilities p* kn {t) and 
Pk,n(t) coincide, so that §5§ holds. ■ 

It should be noted that if the ratio \ n J fi n 
is a constant, a solution of system ([3]) can be 
easily obtained. Indeed, if there exists a pos- 
itive constant c (c ^ 1) such that fi n = c\ n 
for all n € Z, then one has: 



z/ n = l + /?c r ' 
with p > 0. 



(n e Z) 



3. First— passage— time densities 

Let 

= 

= — Pr ( inf{t > : X t = s} < t 1 X = k). 

be the first -passage-time probability density 
function of X t through the state s, condi- 
tioned upon the initial condition Pr{X = 
k} — 1, with k s. In the following The- 
orem we determine the relations existing be- 
tween the first-passage-time densities of the 
birth-and-death processes X t and X t . 

Theorem 2 Under the assumptions of The- 
oremUl the first-passage-time density g k ,s(f) 
of process X t with rates @ is given by 

= — 9k,a(t) {k^s). (6) 



Proof. The transition probabilities and the 
first-passage-time density of X t are related 
by the integral equation 

Pk,n(t) = f gk,s{'&)Ps,n{t-'d)d-&, 

Jo 

which holds as k<s<n or n<s<k. 
Making use of relation ([5]) one obtains: 

Pk,n[t) = — I 9k,s{P)Ps,n(t - ■&) d-d. 

v s Jo 

However, one also has: 

Pk,n{t) = [ g k ,s{fi)Ps,n{t- $)d$. 

Jo 

Comparing the last two integral equations, 
due to the uniqueness of their solution, iden- 
tity ([6]) follows. ■ 

Bearing in mind that a first-passage-time 
density could be "defective", let us now pin- 
point some further relations. If sequence 
{z/ n }„ is increasing, from ([6]) one has: 

• 9k t a(t) > g ktS (t) if s > k; 

• 9k,s(t) < 9k,s(t) if s < k. 

Instead, when {v n } n is decreasing the inequal- 
ities between the first-passage-time densities 
are inverted. 

For all k ^ s let us now introduce the ulti- 
mate crossing probabilities 

/•oo roo 

V k .,= / 9k,s(t) dt i ^k,s= / 9k,s(t) dt - 











Integration of both sides in (jBJ) over (0, oo) 
immediately yields 

Vk,s = —Vk,s (k^s). 

As {u n } n is a strictly monotonic sequence, the 
ultimate crossing probabilities Vk, s and Vk, s 
are never equal; hence, in particular, they can- 
not be simultaneously unity. 



4. An example 

In this Section we give an example of appli- 
cation of previous results. Let X t be the bi- 
lateral birth-and-death process with constant 
rates X n = A > 0, \i n = /i > (n G Z). The 
transition probabilities and the first-passage- 
time density of X t are given by (cf. Abate et 
al P): 



Pk,n(t) = 

and 



I n - k {2t\J\ix) 



t 



(7) 



(8) 



t 



fx 



respectively, where l k is the modified Bessel 
function of first kind: 



h(2z) = £ 



,k+2j 



Note that when A = /i equation does not 
admit of a positive non-constant solution. In- 
stead, if A 7^ ji Theorem CD holds. Further- 
more, from equality ([3]) for all n G Z one has: 



/A, 



l + /?c r ' 



(/3 > 0), 



where c = /x/A. Recalling expressions HI for 
all n G Z the rates of the birth-and-death 
process X t are given by: 



A n — A 



l + pc 



,n+l 



,n-l 



l + /?c' 
1 + /?c n 



Furthermore, making use of equations (JSJ) 
and (j7|), the transition probabilities of V 4 fol- 
low: 



Pk,n(t) 



l + (3c n 
l + pc k 

Xe -(A+ M )t c (fe-n)/2 In k (2tJ)^). 



Finally, making use of Theorem (2], from (jSJ) 
we obtain the first-passage-time density of X t 
from state k to state s: 



9k,s{t) 



l+(3c s \s-k\ 
l + pc k t 

Xe -(A+^)t c (fc-.)/2 Is _ k{ 2tJ\^). 
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